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Abstract 

We study evolution equations for electric and magnetic field amplitudes in a ring cavity 
with plane mirrors. The cavity is filled with a positive or negative refraction index material 
with third order effective electric and magnetic non-linearities. Two coupled non-linear equa- 
tions for the electric and magnetic amplitudes are obtained. We prove that the description 
can be reduced to one Lugiato Lefever equation with generalized coefficients. A stability anal- 
ysis of the homogeneous solution, complemented with numerical integration, shows that any 
combination of the parameters should correspond to one of three characteristic behaviors. 

PACS: 42.65.Sf, 05.45.-a 

1 Introduction 

During the last decade, composite material developments allowed the experimental realization of 
materials with negative refraction index [TJ [2] , theoretically predicted by Veselago [3j . These ma- 
terials have simultaneously a negative dielectric permittivity and negative magnetic permeability, 
a property that is not found in natural materials. Negative refraction index materials (NRM) have 
many new and interesting properties, such as a refracted wave on the same side as the incoming 
wave respect to the surface normal and Poynting vector in the direction opposite to wave-vector. 
A number of applications has been proposed based on their new properties, such as perfect lenses, 
phase compensators, and electrically small antennas [H [5j [6j [7] . 

Most studies on NRM consider linear relations between electric field and polarization, and 
magnetic field and magnetization. The study of nonlinear effects acquired an increasing interest 
during the last years. It is known that the propagation of an electromagnetic wave in a Kerr 
nonlinear material, with positive refraction index, can be described by a nonlinear order parameter 
equation of the Schrodinger type [8] , and that the same kind of equation can be extended to an 
NRM [21 Uni HH H21 US] ■ Nevertheless, it has been shown [TH [TJ] that a composite metamaterial 
with negative refraction index can develop a non-linear macroscopic magnetic response. This means 
that, although the host medium has a negligible magnetic non-linearity, the periodic inclusions of 
the metamaterial produce an effective magnetic non-linear response when the wave-length is much 
larger than the periodicity of the inclusions. In the rest of the paper, when we speak about magnetic 
non-linearity we refer to this macroscopic magnetic response present only in metamaterials. Electric 
and magnetic non-linearities in composite materials have been also analyzed in, for example, [16, 

nzi us us] 

In this paper, we are interested in the analysis of the equations that describe the electric and 
magnetic fields in a ring cavity with plane mirrors containing a material with negative or positive 
refraction index and with electric and magnetic non-linearities. The aim of the work is to obtain 
a simple mathematical description of this system, useful to identify relevant parameters and to 
analyze typical behaviours. 

The paper is organized as follows. In Sect. [5]we present the equations for the evolution of the 
electric and magnetic field amplitudes in a ring cavity (the derivation, starting from two coupled 
non- linear Schrodinger equations, is in the appendix), and show that the description can be reduced 
to one Lugiato Lefever equation 20] with generalized parameters. In Sect. [31 we analyze the effects 
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of dissipation. In Sect. [4] we use a linear stability analysis to identify three typical situations that 
arise depending on the signs of the three parameters of the equation. Numerical integration 
supports and completes the previous analysis. In Sect. [5] we present our conclusions. 

2 Equations in the cavity 

The equations that describe the behavior of the electric and magnetic fields in a plane perpendicular 
to light propagation are of the type of the Lugiato-Lefever (LL) equation. The LL equation is a 
simple mean field model that has been useful for the analysis of pattern formation in a cavity with 
flat mirrors containing a Kerr medium and driven by a coherent plane- wave field, see also [211 122) . 

We will first analyze the problem of free propagation (without mirrors) of an electromagnetic 
wave in the material and afterwards will use the resulting equations to derive the behavior in the 
cavity. 

We will consider a linearly polarized driving field with frequency ujq. Let us suppose that the 
electric field is in the x direction and the magnetic field is in the y direction. The starting point are 
the Maxwell's equations and the constitutive relations for the electric displacement, D = eoE + P, 
and the magnetic induction, B = [IqH + /iqM . (There is an interesting alternative approach, 
described in [23], in which the fields E, D and B are used, with D = eE and B = H, where e is a 
generalized dielectric constant.) 

Considering an isotropic metamaterial with third order non-linearities, the non-linear relation 
between polarization P (in the x direction) and electric field E is 

/oo />oo 
X ^{t-r)E{T)dT + e / X ( eHt-n,t-T 2 ,t-T 3 ) E{ n )E{T 2 )E(T 3 ) dT X dr 2 dr^ (1) 
-oo J —oo 

where Xe and x^ are the linear and nonlinear electric susceptibilities J5J. The same geometric 
arguments used in (QJ can be applied to the relation between magnetization M and magnetic field 
H: 

X$(t-T)H(T)dT+ / x^ ) (t-r 1 ,t-T2 ) t-T 3 )H(n)H(T 2 )H(r 3 )dndT 2 dT 3 , (2) 

-oo J — oo 

with Xrn and Xm being the linear and nonlinear magnetic susceptibilities. Let us remark that Xe , 
Xe 3 "* j Xrn and Xrn describe macroscopic effective response of the metamaterial valid for a driving 
field with a wave-length much larger than the periodicity of the inclusions. We are assuming 
that the field intensities are small enough in order to neglect higher order terms. The effective 
macroscopic magnetic non-linearity can become relevant in metamaterials, as shown in |14j . in 
contrast with what happens in conventional optics, where magnetic non-linearities are negligible. 

The classical multiple-scales perturbation technique will be applied, in which it is assumed that 
the light is quasimonochromatic and can be represented by a plane wave, with frequency loq and 
wave- number fcoi propagating along the z axis and modulated by a slowly varying envelope. The 
envelope depends on space R = (X, Y, Z) and time T variables that have characteristic scales much 
greater than the scales given by 1/fco and 1/cjo- The fields E and H can be written in the following 
way 

E = £(R,T)e l(feoZ -" ot) +c.c. 

H = "H(R,T)e i(fe ° z ~" ot) +c.c. (3) 

where £ and H are the slowly varying amplitudes. Similar relations hold for P and M. 

Let us define the wave-number k(uj) = wn(w)/c; the refraction index is n(ui) = ±-y/e r (w) (i r (uj) 
(it takes the negative sign when both, e r and fi r are negative [3]), where e r (co) — 1 + Xe '(w) and 
fi r (u>) = 1 + M are the relative permittivity and the relative permeability respectively. In 
particular, ko = ui^n/ 'c, with n = n(u>o), and we call k' and k" the derivatives of fc(w) evaluated in 

LU = UJq. 
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Figure 1: Ring cavity with the nonlinear optic material of size L. 



The details of the application of the multiple-scales technique in our case are essentially the 
same to the ones described in Ref. [8], Sect. 2k, for a positive refraction index material with only 
electric non-linearity. After this process, we arrive to the following coupled non linear Schrodinger 
equations for the envelopes of the electric and magnetic fields: 
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where, in the left hand side, the transformation £ = z, r = t — k'z, was used; the transverse 
Laplacian Vj_ = {~§£Si~§gz) is defined in the plane perpendicular to the z axis. The relative 
permittivity and the relative permeability, e r and /v, are evaluated in uiq. The non-linear electric 
and magnetic susceptibilities, Xe and Xm , are the Fourier transforms evaluated in (u>o, uiq, — ujq). 
Dissipation (linear or non-linear) is neglected, so that x& \ Xm , and Xm are real quantities. 
This is an often met approximation in conventional optics, but dissipation could play a relevant role 
in metamaterials. In this section we will derive the equations for a material without dissipation, 
and in Sect. [3] we will analyze how these equations are modified when dissipation is taken into 
account. 

The nonlinear electric susceptibility is usually written as Xe 3 '' = ot/E 2 , where a = ±1 stands 
for a focusing or defocusing nonlinearity and E c is a characteristic electric field. Eqs. (@| and (JS|) 
are an extension of the analysis performed in [SJ [TD] to include the magnetic non-linearity, and are 
similar to the ones derived in [13] . 

For an NRM, fco is negative, but this does not modify the sign of the non-linear terms in Eqs. 
|4|) and since ko/e r and fco/Mr are always positive. The difference appears in the diffraction 
term, that becomes negative for an NRM. 

Fig. [T] shows the geometry of the ring cavity. At the entrance mirror, the input field has 
amplitudes E ln and H m , with H m = E- ln yeo/Mo- The non-linear material has length L, the 
transmission coefficient in the right end of the material for the electric field amplitude is t+ = 1 + C 
and for the magnetic field amplitude is t_ = 1 — £, where £ = (rj — ?7o)/( 7 ? + Vo) with rjo = \/ e o/ Mo 
and r\ = y/e/(i. The transmission and reflection coefficients in the input mirror are fj and r, . The 
time of a round trip is T r , the phase accumulated by the wave in a round trip is <f>, and the detuning 
between the pump field and the cavity mode is 9 = <f> mod 27r, with <C 1 since we consider that 
the cavity is close to resonance. It is convenient to define the quantity p = |r,|(l — £ 2 ), the length 
I = y/Lp/[2\ko\(l — p)], and the characteristic electric field C — t + E c ^/2e r (l — p)j (3Lk p). Using 
the non-linear Schrodinger equations (UJ) and ([5]) for the fields inside the cavity, and after the 
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following change of variables: 

Ai = —, A 2 = —*—, x =x/l, y =y/l, t =t—= — , (6) 

we arrive to two Lugiato-Lefever type equations (the details of the derivation can be found in the 
appendix) : 

fl A 

' 1 An - (1 + ie)A, + i^\A x + l (a\A 1 \ 2 + 7 |A| 2 )A (7) 



Of 

" 2 - A in -(l + iO)A 2+i pviA 2 + i(a\A 1 \ 2 + 1 \A 2 \ 2 )A 2 (8) 



dt 

where = —9p/(l — p), (3 = ±1 is the sign of the refraction index n, 7 = Xm l\Xe^ I; an d 
An = 7j^y-Ei n /C. The primes have been omitted to simplify the notation. 

Let us analyze the difference of the fields: R = A — A 2 . It is possible to prove that R decays 
exponentially to zero. The equation for R is 

dR 

— = -(l + iO)R + i(3VlR + i( a \A 1 \ 2 +-f\A 2 \ 2 )R, (9) 

from which we get 

= ~2\R\ 2 + ipRV 2 L R-if3RV\R, (10) 

ot 

where R is the complex conjugate of R. We use the expression of R in terms of its Fourier transform, 
R{x,y,t) = -^Kp; J dk x dkyR kx ^ ky e i( ' k * :X+kyV \ The last two terms in Eq. <jTU]) become 

RV\R- RV 2 X R = 

J dk x dk y dq x dq y R k ^ ky K, qy Z l[k *- q * )X e^-^iql + q 2 y - k 2 x - k 2 y ). (II) 

Let us consider the transverse average of Eq. (jXTJ) , defined as the integral over the plane x-y, 
(RV 2 j_R - RV 2 ± R) = 

J dk x dk y dq x dq y R k ^ ky R q ^ v (q 2 x +q 2 y -k 2 x -k 2 y ) J dx dy e^~^ 

= J dkx dk v d 1 x d 1y R k*,k y Rq x , qy {ql +q 2 y ~k 2 x - kl) 8{k x - q x ) 8{k y - q y ) 

= 0. (12) 



Therefore, from Eq. (| 10[) we get 

^ = -2(|i?| 2 ), (13) 

that means that the transverse average (|-R| 2 ) decays to zero exponentially with a characteristic 
time 1/2 (or T r /[2(1 — p)] for the previous time scale). If, after a transient, we have (|-R| 2 ) = 0, 
then, since |i?| 2 > in every point of the transverse plane, we have that R = in every point. 
Then, after a time of order 1/2, we have that A\ = A 2 . In this situation, the description of Eqs. 
([7]) and © is further simplified to 

8 A 

— = An - (1 + iQ)A + ipV 2 ± A + i~f'\A\ 2 A, (14) 
ot 

where A — A = A an d 7' = a + 7. Eq. (TT4")) has the same form than the Lugiato Lefever equation, 
but there are some differences. In the original version, that considers only an electric non-linearity, 
it was shown |20j that the sign of the detuning must be equal to the sign of the non-linear 
coefficient, 7', but this is not necessarily the case in Eq. pT|) . In addition, in (Tl"4"f the sign of the 
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diffraction term, /3, can be negative or positive depending on the material being an NRM or not, 
and the magnetic non-linearity is included in the coefficient 7'. 

In the particular case in which there is only an electric non-linearity, i.e. 7 = 0, and using 
the result of [20] that says that the sign of Q is equal to a, that is equal to 1 (—1) for a focusing 
(defocusing) electric non-linearity, we get, 

— = A in - (1 + ia\G\)A + i/3V\A + ia\A\ 2 A. (15) 

Eq. (|15[) presents an interesting symmetry. Taking its complex conjugate, and assuming that A[ n 
is real, it can be seen that a focusing non-linearity and a positive refraction index material (a = 1 
and f3 = 1) is equivalent to the defocusing and NRM case (a — —1 and (3 = —1). Also the case 
(a = 1,/J = —1) is equivalent to (a — —1,(3 — 1). As we will see below, the equivalences are more 
involved when the magnetic non-linearity is considered, since the signs of the three coefficients (0, 
/3 and 7') are in general not related among them. 



3 Dissipation 

Although it is common to neglect dissipation in conventional optics, this is not in general an appro- 
priate approximation for a metamaterial. In this section we analyze how the previous equations are 
modified when dissipation is taken into account. In this case, the electric permittivity and magnetic 
permeability are complex quantities. The wave- number that was defined as k(u) = t r (ui)p r {ui) 
will also have an imaginary part. At the frequency ujq we have 

fc(w ) = fc + ikj (16) 

where fco is the wave- number of the plane wave modulated by the slowly varying amplitude ([3]), 
and ki is the imaginary part. 

We will consider that the attenuation distance is much larger than the wave-length. This 
means that ki -C ko. In the Gigahertz range, it is possible to build an NRM with small (and even 
negligible) imaginary parts of e r and fi r as was shown in, for example, [2H [25] . 

It is assumed that fcj is of second order in the small parameter used in the multiple-scales per- 
turbation technique, see [8l p. 101]. The result is an additional term in the non-linear Schrodinger 
equations for the envelopes of the electric and magnetic fields 

f - 4^ivi* + ^(|W + f *,.)«-**. « 18 , 

The additional terms (—kjE and —kjTi) represent an exponential decay of the amplitudes due to 

_____ /g-j /g\ 

dissipation. In Eqs. (fTT|) and (fI5]l. the non-linear coefficients \e /e r and Xm / can be taken as 
real quantities since the contribution of the imaginary parts correspond to terms of higher order 
in the expansion of the small parameter. 

Following the steps indicated in the appendix, we arrive to a couple of Lugiato-Lefever equations 
that have the same form that Eqs. and (jl]), where now the change of variables is given by 

X ' = X /l>, y> = y/l', t , = t (l-P + PLk I ) > (19) 




with C = t + E cA /2e r {\ - p + pLkj) / '(3Lk Q p) and I' = y/Lp/[2\k \(l - p + pLkTj\. The scaled 
detuning is now defined as = —Bp/ (1 — p + pLkj), and the input field is A- m = (!_p_[_p^fc f ) E in /C. 
The definitions of the rest of the coefficients that appear in Eqs. (JTJ) and ([5]) remain the same. 

Therefore, the case of small dissipation can still be correctly described by Eqs. (JT]) and {5}. 
It is not difficult to see that the demonstration of the previous section that the amplitudes A\ 
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and A2 are equal after a transient still holds for this dissipative case. Consequently, the simplified 
description of Eq. (fl4|) also holds. 

The situation is different if the attenuation distance is of order of the wave-length, i.e., ki ~ fco, 
that corresponds to a factor of merit of order one. This is generally the case for optical frequencies. 
Now, the assumption that light can be represented by a plane wave modulated by a slowly varying 
amplitude, Eq. ([3]), is no longer valid, because the decay of the fields due to dissipation takes place 
in a fast space scale, and this decay can not be described by the envelopes £ and H. The situation 
when the length of the material, L, is much greater than the wave-length, is not interesting since no 
field will be detected at the output. There are many proposals to reduce losses in metamaterials, 
some of them are [2i [571 HH1 [23 1201 E] • 

4 Stability of homogeneous stationary solutions 

The homogeneous stationary solutions, Aq, of Eq. (Tl4|) are obtained by solving 

A in = (l+ie-iYl )A , (20) 

where Iq = \Aq\ 2 . For the intensities we have 

7 in = (1 + e 2 )/ - 20 7 7 2 + 7 ,2 / 3 , (21) 

where /;„ = |A; n | 2 is the intensity of the input field. It is well known that Eq. ([21]) presents 
bistability for |8| > -s/3, i.e., there is a range of values of /;„ for which there are two stable 
solutions of Iq. Note that, since = —9p/(l — p), and 9 <C 1, bistability can only be attained 
when p = \ri\(l — ^ 2 ) ~ 1, i.e., when the entrance mirror has a high reflectivity and the material 
has a high transmissivity. Another condition to have bistability is ©7' > 0, this condition is 
automatically fulfilled when there is only an electric non-linearity, i.e., when 7 = and = a|0|. 
In terms of the solutions for the intensity Iq obtained from (j2TJ) , the homogeneous solution is 

A o = (l-i0 + i 7 7o)/o/A n , (22) 

where we have assumed, without loss of generality, that A ln is real. 
A linear stability analysis of Aq gives the following eigenvalues, 

A± = -1 ± #0 - /3fc 2 - 0][/3fc 2 + - 7 ' Jo], (23) 

where k is the wave-number of a small perturbation. From (|23p we can draw the marginal stability 
curves of Fig. [2J The same figure applies to both cases: negative and positive refraction index, 
that correspond to the negative and positive parts, respectively, of the horizontal axis. In the 
regions enclosed by the curves, the homogeneous solution Ao becomes unstable. Both regions do 
not appear simultaneously: one corresponds to 7' > and the other to 7' < 0. A value of = 
was used in the figure, but a different value of only represents a horizontal shift of the curves. 
The homogeneous solution becomes unstable when the intensity is increased from zero an reaches 
the value Iq | 7 '| = 1. Over the instability threshold, and close to it, it is known that a hexagonal 
pattern appears (32] [33], when the critical wave number is different from (this happens when 
there is not bistability or, more strictly, when /30 < 2). 

Eq. (|14j) and its complex conjugate represent the same physical situation. The only difference 
present in the complex conjugate equation is the sign in front of the three real parameters 0, /3 
and 7'. Using this equivalence we can identify three typical situations: (a) /3 7 ' > and any value 
of 0; (b) ^ < and 7 ' < 0; and (c) /3 7 ' < and 7 ' > 0. 

In case (a) /3 and 7' have the same sign. We can see in Fig.[2]that, for j3 = — 1 or 1, the marginal 
stability curve is present, so that, as the input field intensity is increased, the homogeneous solution 
becomes always unstable. The critical values of field intensity and wave-number are given by, 

7 0)C | 7 '| = 1 and fc 2 = 2 - |0| for /30 < 2 
io, c |7'| = 2|0|/3 - V© 2 - 3/3 and k c = for (3Q > 2 
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Figure 2: Marginal stability curves, Iq | 7 '| versus (5k 2 , for 8 = 0. In the dashed region the 
homogeneous solution becomes unstable. 

In case (b), the homogeneous solution is always stable. Bistability is not possible in this case 
since 87' < 0. 

In case (c), the homogeneous solution becomes unstable only if there is bistability, i.e., when 
|8| > -\/3- The critical values are, 



If 8 > y/3, and 7' > 0, the vertical axis in Fig. [2] is shifted to the right and crosses the 
marginal stability curve in two points that give a range of values of Io for which the solution Ao 
becomes unstable under homogeneous perturbations. This corresponds to the unstable branch 
when bistability is present. Fig. [3] shows a sequence of plots with the marginal stability curve 
with the values of the homogeneous solutions indicated on the vertical axis for increasing values of 
the input intensity I m and for 8 = 3 and 7' = 1. As I- m increases from 2.9 to 2.95 a saddle node 
bifurcation takes place in the upper branch of the marginal stability curve (the bifurcation happens 
in the two dimensional phase space composed by the real and imaginary parts of Ao). As I[ n is 
further increased, the value of the intensity of the homogeneous solution of the unstable branch 
(the one that is inside the instability region delimited by the curve) decreases, until it merges, in 
another saddle node bifurcation, with the lower homogeneous solution. 

We have an example of case (c) in Fig. [3] for (3 = —1, where the lower homogeneous solution 
becomes unstable before it merges with the homogeneous solution of the unstable branch (see 
dotted line in bottom right plot of Fig. [3J. The saddle node bifurcation in the lower branch takes 
place at Iq approximately equal to 1.18 for the parameters of Fig.[3j but a modulational instability 
with a wave- number different from zero takes place at a lower value of Iq (Jq = 1). 

Fig. |4] shows four snapshots of numerical integration results for this case. The initial condition 
is the lower homogeneous solution slightly above the instability threshold (Ao = 0.457 — iO. 908, i.e. 
Io = 1.03), plus noise of amplitude 0.005. Eq. (fl"4"|) was integrated in a square domain of 256x256 



I jC | 7 '| = 2|8|/3 - \/<d 2 - 3/3 and k c = for V3 < |8| < 2 
Z 0iC | 7 '| = l and kl = \Q\-2 for |8| > 2 
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Figure 3: Marginal stability curve, Iq versus f3k 2 , for the bistability case with 6 = 3, and 7' = 1. 
The sequence corresponds to increasing values of ii n , from left to right and from top to bottom: 
L m = 2.9, 2.95, 4 and 5.05. The small circles on the vertical axis show the positions of the 
homogeneous solutions. The negative horizontal axis corresponds to case (c) and the positive one 
to case (a). 



points with periodic boundary conditions, using a Fourier series representation and a 4th order 
Runge-Kutta temporal scheme for the non-linear terms. At time t — 180 the field still appears 
homogeneous. At t = 270 a pattern characterized by the critical wave number appears (k c — 1). At 
t = 297 we can see some spatial domains where the field takes the value of the upper homogeneous 
solution. These domains grow until the whole system becomes homogeneous. Therefore, the 
numerical results indicate that, as the input intensity is increased, the lower homogeneous solution 
becomes unstable, but this instability does not give rise to a stable non-homogeneous pattern. 
Instead, the system evolves to the upper homogeneous solution, that is always stable for /3 = — 1. 

Continuing with the analysis of Fig. [31 the situation is different for j3 = 1, that corresponds 
to case (a); in this case the lower homogeneous solution is always stable until the saddle node 
bifurcation takes place. Therefore, as the input intensity is increased, the lower homogeneous 
solution does not simply cross the marginal stability curve and lose its stability, but, instead, 
it does no longer exist as a stationary solution since the homogeneous mode, k = 0, becomes 
unstable. We performed numerical integration of Eq. (1141) for an input intensity slightly above 
the value corresponding to the saddle node bifurcation. For different values of the homogeneous 
initial condition, the evolution is qualitatively similar to the one shown in Fig. Initially, only the 
modes close to are unstable, so that, for short times, the value of the intensity increases keeping, 
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Figure 4: Four configurations of \A\ 2 against position in the transverse plane, obtained from 
integration of Eq. (fT4")l for increasing values of time. Parameters are /3 = — 1, 9 = 3 and 7' = 1, 
corresponding to case (c). The initial condition has intensity Iq = 1.03 (/;„ = 5.02), slightly above 
the critical value (Iq. c = 1). The system evolves from the lower to the upper homogeneous solution. 
The space and time integration steps are Ax = 0.25 and At = 0.003. The grey scale is logarithmic, 
black corresponds to 0.4 and white to 4.5. 

approximately, the homogeneous shape of the field. For larger times, the system evolves to a state 
of optical turbulence or spatiotemporal chaos. 

As mentioned before, in case (a) without bistability, an hexagonal pattern appears close and 
above the instability threshold. It was shown in Ref. t 34 that, as the input intensity is increased, 
there is a sequence of different spatiotemporal regimes: oscillating hexagons, quasiperiodicity, 
temporal chaos and optical turbulence. When there is bistability, the state of optical turbulence 
is reached directly as a transition from the homogeneous state as the input intensity is increased 
(see Fig. [5]). 

The differences between the stability of homogeneous solutions for positive (j3 — 1) or negative 
(f3 = —1) refraction index materials, when bistability is present, is clearly presented in Fig.|6l The 
figure shows the curve Iq against /; n . The shape of the curve is the same for /3 = 1 or /3 = — 1, but 
the stability ranges are different. 

The stability analysis is similar to the one presented in [9] , where only an electric non-linearity 
was considered. The inclusion of the magnetic non-linearity is explicitly represented by factor 7 in 
(|14p . In addition, a result that, to our knowledge, was not previously reported, is the equivalence 
between positive and negative refraction index materials when the signs of the detuning and 
the non-linear coefficient 7' are changed. This symmetry is illustrated, for example, in Fig. [2J for 
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Figure 5: Intensity \A\ 2 for increasing values of time calculated from Eq. (Q~4|). Parameters are 
/3 = 1, 8 = 3 and 7' = 1, corresponding to case (a). The homogeneous initial condition has 
intensity equal to 1.19 plus noise of amplitude 0.005, and the input field has /;„ = 5.11. The final 
state corresponds to spatiotemporal chaos that is reached independently of the value of the initial 
homogeneous state. Space and time integration steps: Ax — 0.125 and At = 0.001. The gray scale 
is logarithmic with black equal to 0.1 and white equal to 60. We used a different scale for t = 22: 
black=3, whitc=4, in order to show the transient emerging pattern. 
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Figure 6: Intensity of the homogeneous solution Iq against intensity of the input field I ln in the 
bistability case, with = 3 and 7' = 1. The continuous curve corresponds to an always stable 
homogeneous solution. In the rest of the curve, the stability range depends on the sign of the 
refraction index, j3. The part with negative slope (dotted curve) is always unstable for j3 = 1 or 
— 1. The upper part (dash-dotted curve) is unstable only for /3 = 1 but it is stable for an NRM 
(f3 — —1). Equivalently, for the region 1 < Iq < 1.18 the homogeneous solution is unstable for 
/3 = — 1 but is stable for (3 = 1. 



6 = 0. In this case, the behavior for j3 = 1 and 7' > is equivalent to the behavior for j3 = — 1 
and 7' < 0. 



5 Conclusions 

Starting from the Maxwell equations, with non-linear polarization and magnetization, it is possible 
to obtain, using the multiple-scales technique, two coupled non-linear Schrodinger equations for 
the electric and magnetic field amplitudes. From these equations we derived the evolution of the 
fields in a ring cavity with plane mirrors containing a material with positive or negative refraction 
index and with effective electric and magnetic non-linearities. We proved that the description can 
be reduced to only one equation that has the same form of the Lugiato Lefever equation [20] for a 
positive refraction index material with only electric non-linearity. 

An original contribution of the paper is the generalization of the Lugiato Lefever equation to 
electric and magnetic non-linearities. Let us note that, as was shown in [14], effective macroscopic 
magnetic non-linearities can become relevant in composite materials that are used to generate 
a negative refraction index. The resulting equation has three real parameters: the sign of the 
refraction index /?, the detuning 0, and the non-linear coefficient 7'. In the original version of the 
Lugiato Lefever equation, the sign of the detuning is equal to the sign of the non-linear coefficient, 
corresponding to the self-focusing or self-defocusing cases. In the present version, both signs are 
independent; in addition, the diffraction coefficient can be positive or negative depending on /?. 
This generalization of the parameters allows the equation to represent new physical situations. 
Despite the fact that more free parameters necessarily makes the analysis more complex, using a 
linear stability analysis we have shown that any combination of the parameters must correspond 
to one of only three typical cases. In case (b), ^7' < and 87' < 0, the homogeneous solution is 
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always stable. In case (c), {5~f < and 87' > 0, the homogeneous solution can become unstable 
only when there is bistability, but, in this case, numerical integration shows that the final state is 
the homogeneous solution of the upper branch. Only in case (a), ^7' > 0, the destabilization of 
the homogeneous solution, as the input field intensity is increased, gives rise to a non-homogeneous 
state. In this last case, if there is not bistability it is known that, close to the instability threshold, 
the asymptotic state is an hexagonal pattern. If there is bistability, numerical results for = 3 
show that there is a transition from the homogeneous state to optical turbulence as the input field 
is increased. 

In summary, the description, in the plane perpendicular to propagation, of the evolution of 
the electromagnetic field in a cavity with a non-linear material with positive or negative refraction 
index has been reduced to a Lugiato Lefever equation with three parameters: ©, 7' and /3. These 
quantities are functions of the much larger set of parameters of the original description (based on 
the Maxwell equations). One of the aims of the work was the identification of relevant parameters 
since they allow a better understanding of typical behaviours that the system can develop. This 
kind of analysis is much more difficult to perform with the original description. 
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Appendix 

In this appendix we present the derivation of the cavity equations ([7]) and ((8|) starting from the non 
linear Schrodinger equations (HJ) and ([5]). It is essentially an extension of the procedure described 
in [9] to the case of electric and magnetic non-linearity. 

Using Eqs. (j3| and (JSJ, for a small distance we can write 

T± (£ + *0 = T± (0 + 1 8t NF ± (0 ~ e 4 5 « * T± (£) (24) 
where J r + = £, T- = H, and N is the operator 

k" B 2 1 1k„ /v (3) v (3) \ 

N = -TW + ^ + -T[^ 2 + X t m2 ) (25) 

Fig. [T] shows the scheme of the cavity with the non linear material of length L whose left and 
right ends are at positions a and b. We call a + and a_ the positions immediately to the right 
and to the left of a respectively, and similarly for b + and 6_ . Using the impedance of free space, 
Vo = \Ao/Vo, an d of the material, 77 = v/i//i, the transmission coefficient for the electric field 
amplitude, when light goes from 6_ to 6+, is 

t+ = = 1 + C (26) 

where £ = (77 — T]o)/(r] + T}q) is a small quantity, so that the transmission coefficient is close to 
1. (Even when 77 is not close to 770, a transmission coefficient close to 1 can be achieved by filling 
the cavity with a substance with an impedance close to 77.) The transmission coefficient for the 
magnetic field amplitude in b is, 

t= Jv^ = 1 _ c (2?) 



12 



Using also the transmission and reflection coefficients in the input mirror, t{ and rj, we obtain the 
following relations 

F ± (b+) = *±J±(6_) 

7±(a_) = |ri|e^.F±(& + )+*i^n± (28) 

where F- m+ = E- m , F- m - = H- ln , and 4>, the phase accumulated in a round trip, includes the phase 
change of 4ir due to reflection in the four mirrors. We will consider that the system is close to 
resonance so that the detuning is 9 <C 1, with 9 = <fi mod 2ir. Let JF*± be the value of T± after 
one round trip in the cavity, that takes a time T r . Applying the relations (|24p and (f2"5f we find 
(subindices ± are removed for simplicity) 

F« = UF m + \ n \(l - Q 2 y 6 e lL *F ~ tiFin + \n\(l - C 2 )(l + + iLN)F. (29) 

According to (|23|) . the operator iV depends on the fields inside the material, but, in Eq. (|2T)1) . J 7 
corresponds to the fields outside the material, so N should be evaluated in t T F±. Writing the 
time derivative of F as (F* — F)/T r , and defining p = |rj|(l — £ 2 ), we obtain, 

Tr ~dt = kFin + (P ~ l) T + lpBT + NT - (30) 

After replacing the operator N by its definition (|25p , and using the change of variables defined by 
Eqs. ([6]), we arrive to Eqs. (J7]) and ([SJ for the amplitudes of the electric and magnetic fields in the 
cavity. It can be shown that, in the cavity equations, the dispersion term proportional to k" in 
(|23|) can be neglected. 

References 

[1] D.R. Smith, W.J. Padilla, D.C. Vier, S.C. Nemat-Nasser and S. Schultz, Phys. Rev. Lett. 84 
4184 (2000) 

[2] R.A. Shelby, D.R. Smith and S. Schultz, Science 292 77 (2001). 

[3] V.G. Veselago, Usp. Fiz. Nauk 92 517 (1967); V.G. Veselago, Sov. Phys.-Usp. 10 509 (1968). 
[4] Z. Jaksic, N. Dalarsson and M. Maksimovic, Microwave Review 12(1), 36-49 (2006). 
[5] S.A. Ramakrishna, Rep. Prog. Phys. 68, 449 (2005). 

[6] A. Grbic and G.V. Eleftheriades, Physical Review Letters 92, 117403 (2004). 
[7] N. Fang, H. Lee, C. Sun and X. Zhang, Science 308, 534-537 (2005). 

[8] J.V. Moloney and A.C. Newell, Nonlinear Optics (Westview Press, Boulder, Colorado, 2003). 

[9] P. Kockaert, P. Tassin, G. Van der Sande, I. Veretennicoff and M. Tlidi, Phys. Rev. A 74, 
033822 (2006); P. Tassin, L. Gelens, J. Danckaert, I. Veretennicoff, P. Kockaert and M. Tlidi, 
Chaos 17, 037116 (2007). 

[10] P. Tassin, G. Van der Sande, I. Veretennicoff, M. Tlidi and P. Kockaert, Proc. SPIE 5955, 
59550X (2005). 

[11] S. Wen, Y. Wang, W. Su, Y. Xiang. X. Fu and D. Fan, Phys. Rev. E 73, 036617 (2006). 

[12] M. Scalora, M.S. Syrchin, N. Akozbek, E.Y. Poliakov, G. D'Aguanno, N. Mattiucci, M.J. 
Bloemer and A.M. Zheltikov, Phys. Rev. Lett. 95, 013902 (2005). 



13 



[13] N. Lazarides and G.P. Tsironis, Phys. Rev. E 71, 036614 (2005); I. Kourakis and P.K. Shukla, 
Phys. Rev. E 72, 016626 (2005). 

[14] A. A. Zharov, I.V. Shadrivov and Y.S. Kivshar, Phys. Rev. Lett 91, 037401 (2003). 

[15] M. Lapine, M. Gorkunov, and K. H. Ringhofer, Phys. Rev. E 67, 065601 (2003). 

[16] S. OBricn, D. McPeakc, S. A. Ramakrishna, and J. B. Pendry, Phys. Rev. B 69, 241101 
(2004). 

[17] I.V. Shadrivov and Y.S. Kivshar, J. Opt. A: Pure Appl. Opt. 7 68-72 (2005). 

[18] I.V. Shadrivov, A.B. Kozyrev, D.W. van der Weide, and Y.S. Kivshar, Optics Express 16 
20266 (2008). 

[19] M. Lapine and M. Gorkunov, Phys. Rev. E 70, 066601 (2004). 

[20] L.A. Lugiato and R. Lefever, Phys. Rev. Lett. 58, 2209 (1987); L.A. Lugiato, W. Kaige and 
N.B. Abraham, Phys. Rev. A 49, 2049 (1994). 

[21] W.J. Firth, A.J. Scroggie, G.S. McDonald and L. Lugiato, Phys. Rev. A 46, R3609 (1992); 
J.B. Geddes, J.V. Moloney, E.M. Wright and W.J. Firth, Opt. Commun. Ill, 623 (1994). 

[22] M. Hoyuelos, P. Colct, M. San Miguel and D. Walgracf, Phys. Rev. E 58, 2992 (1998). 

[23] V. M. Agranovich, Y. R. Shen, R. H. Baughman, and A. A. Zakhidov, Phys. Rev. B 69, 
165112 (2004). 

[24] R. Liu, A. Dcgiron, J. J. Mock and D. R. Smith, Appl. Phys. Lett. 90, 263504 (2007). 

[25] Wang Jia-Fu, Qu Shao-Bo, Xu Zhuo, Zhang Jie-Qiu, Ma Hua, Yang Yi-Ming and Gu Chao, 
Chinese Phys. Lett. 26, 084103 (2009). 

[26] G. Dolling, M. Wegener, C. M. Soukoulis, and S. Linden, Opt. Express 15, 11536-11541 (2007). 

[27] V. M. Shalaev, T. A. Klar, V. P. Drachev, and A. V. Kildishev, Optical Negative-Index 
Metamaterials: From Low to No Loss, in Photonic Metamaterials: From Random to Periodic, 
Technical Dgest (CD) (Optical Society of America, 2006), paper TuC3. 

[28] A. K. Popov and V. M. Shalaev, Optics Letters 31, 2169-2171 (2006). 

[29] P. Tassin, Lei Zhang, Th. Koschny, E. N. Economou, and C. M. Soukoulis, Phys. Rev. Lett. 
102, 053901 (2009). 

[30] J.A. Schullcr, R. Zia, T. Taubner and M.L. Brongersma, Phys. Rev. Lett. 99, 107401 (2007). 

[31] Zhu Wei-Ren and Zhao Xiao-Peng, Chin. Phys. Lett. 26(7), 074212 (2009). 

[32] A.J. Scroggie et al, Chaos, Solitons & Fractals 4, 1323 (1994). 

[33] M. Tlidi, R. Lefever, and P. Mandcl, Quantum Semiclassic. Opt. 8, 931 (1996). 

[34] D. Gomila and P. Colet, Phys. Rev. A 68, 011801 (2003). 



14 



